fig4_data = readtable('fig4.csv');

figure('position',[180,135,450,250]);

plot(fig4_data{:,1},fig4_data{:,2},...
    'o','MarkerSize',4.5,'linew',1,'Color',[0 0.447 0.741],...
    'MarkerEdgeColor',[0 0.447 0.741],'MarkerFaceColor',[0.0745,0.6235,1]);
set(gca,'fontsize',13,'position',[0.1611    0.2492    0.7750    0.6758],'xtick',-60:20:60);
hold all;
ylim(1+0.23*[-1,1]);
xlim(50*[-1,1]);
xlabel('Larmor frequency \omega_a [Hz]');
ylabel('$|1-\mathcal{C}|$','interpreter','latex');
grid on;

plot(fig4_data{:,3},fig4_data{:,4},'-','linew',1.2,'color',[126,47,142]/255);
plot(fig4_data{:,5},fig4_data{:,6},':','linew',2,'color',zeros(1,3));

legend('experiment','simplified model','numerical model');
